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Abstract 

We consider the quantum sinh-Gordon model in this paper. Form fac- 
tors in this model were calculated in gf-|]. We sum up all contribu- 
tions of form factors and obtain a closed expression for a correlation 
function. This expression is a determinant of an integral operator. 
Similar determinant representations were proven to be useful not only 
in the theory of correlation functions |jT^] , but also in the matrix 
models 11-0]. 
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1 Introduction 



The theory of massive, relativistic, integrable models is an important part of 
modern quantum field theory |l7|-p5|. Scattering matrices in these models 



factorize into a product of two-body S-matrices |17j]. Form factors can be 
calculated on the basis of a bootstrap approach |17|-p5 . 



The purpose of this paper is to calculate correlation functions. As usual 
correlation function can be represented as an infinite series of form factors 
contributions. In this paper we sum up all these contributions and obtain a 
closed expression for correlation functions of local operators ( |4.8|) . The idea 
of this summation is the following. We introduce an auxiliary Fock space 
and auxiliary Bose fields (we shall call them dual fields). These fields help 
us to represent the form factor decomposition of a correlation function in a 



form similar to "free fermionic" case. This approach was developed in [26 



[ [27j] and ||. Finally a correlation function is represented as a vacuum mean 
value (in the auxiliary Fock space) of a determinant of an integral operator 
This representation was proven to be useful ||, [§], [|TTJ. It helps 



for asymptotical analysis of quantum correlation functions. Among other 
things this approach helped to calculate asymptotic of time and temperature 
dependent correlation function in Nonlinear Schrodinger equation [^J. 

In this paper we consider the sinh-Gordon model. It is the model of one 
(real) relativistic Bose field <fi in 2D. The action is 



S 



d x 



1 Tfi 

-{d^(x)f §-coshg(j)(x) 

2 g 2 



XI) 



with Z 2 sym- 



It is the simplest example of the affine Toda field theories 
metry — > —<p. The model has only one massive particle. The two-body 



scattering matrix [IS], [29| is given by an expression 



tanh £ 



inB 
2 



tanh U/3 + ^ 



1.2) 



B 



2g 2 



87r + g 2 ' 

We shall consider a real g, which corresponds to a positive B. Later we shall 
use a variable 

x = e 13 (1.3) 
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instead of rapidity (3 . 

We shall use a representation for form factors found in fij-[||] (another 
representation can be found in El) 



!)•••> Hn 



fin 



(0\O(0)\f3 h ... 7 f3 n ) = H n Q n (s)f[ 



i>j Xi + x j 

Here = fii — (3j. A function F min (/3) is holomorphic for real (3 

F min {f3)=N{B)E{i3), 

where 

dx sinh (jxB^J sinh (j;x(l — sinh \x 



;i.4) 



;i.5) 



l(fi) = exp 



AA(5) = exp 



o x 



sinh a; 

d x sinh (j2lE>) sinh (|x(l — ^-£>)) sinh | 



2- .2 x fi 
sin — 
\2tt. 



x 



X 



sinh x 



1.6) 
1.7) 



and (3 = in — (3. The function F min (/?) has a simple zero at (3 = and no 
poles at the strip < Im/? < ir. At /3 — > oo it goes to one: F min (/3) — > 1. 
Functions Q„(xi, . . . , x n ) are symmetric polynomials of variables Xi, . . . ,x n 
given by fl 



Q n (s) = det n _iMy(s y , 



i,j = l,...,n- 1, 



;i.s) 



where 



M J -(s) = <T ai - J -[z-j + s]. (1.9) 

Let us explain the notations. Here and later we suppress the dependency of 
Q n (s) on variables Xj. The index n — 1 in the expression det n _i denotes the 
dimension of the matrix Mij(s). The functions <Jk are elementary symmetric 
polynomials of the fc-th order of the variables Xi, . . . , x n : 



<j k = a k '(x 1 , 



•^n) / j ■^11^12 ■ ■ ■ ^i^i 

h<i2<---<ik 



;i.io) 



and at = if k < or k > n. Here we also suppressed the dependency of 
on Xj. The symbol [m] is a "g- number" defined by 

\n>\ —^ = q - (1.11) 



sin ■ 

g 

sinf 



q-q- 
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where q = exp{iB /2}. The number s in ( |1.8| ), ( |1.9| ) is an arbitrary integer, 



depending on specific choice of an operator O in (|T1 

Finally, the constants in (|1.4|) are normalization constants 



4 s i n 

H2n+i = Hif/ 1 , H 2n = Hof/ 1 , /i = — — t— r, (1-12) 

-*min \1Tt ) 

where H and Hi also depend on specific operator O. For instance, the form 



factor of local field is given by (L4) with s = and 

#o = (O|0(O)|O) = O, 

1 (1-13) 
Hi = (O|0(O)|/5) = - 7 =. 

A correlation function of an operator O can be presented as an infinite 
series of form factors contributions 

oo . Jn Q n 

<0|O(0,0)O(x,0|0> = E/ ZJ7^\ F n^---^n)\ 2 \{e~ mrco ^. (1.14) 

In the present paper we sum up this series explicitly. Now let us present a plan 
of the paper. Section |2| is devoted to a transformation of the determinants 
(O), ( |1.9| ) to a form, which is convenient for summation. In section 3 we 
introduce auxiliary quantum operators — dual fields — in order to factorize an 
expression for a correlation function and to represent it in a form similar to 



"free fermionic case" . In section 4 we sum up the series ( 1.14j) into a Fredholm 



determinant. In section 5 we use the Fredholm determinant representation 
for derivation of an asymptotic behavior of correlation functions. 

2 A transformation of the form factor 

A determinant of a linear integral operator I + V can be written as 



det(I + V) = T [ n det 



/ V(xi,xi) ■ ■ ■ V(x u x n ) \ 
dxi ■ ■ ■ dx n , , V(x 2 ,x 1 ) ■■■ V(x 2 ,x n ) 

V V{x n ,xi) ■ ■ ■ V(x n ,x n ) J 



■ (2-1) 
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Thus, in order to obtain a determinant representation for correlation func- 
tions one need to represent the form factor expansion ( 1.14j) in the form 



( |2.1p . Determinants of integral operators, which we consider also can be 
called Fredholm determinants. 



The form factors (|1.4j ) are proportional to polynomials Q n (s), which in 



turn are equal to determinants of (n — 1) x (n — 1) matrices (|L 

Q n (s) = deWM^s). (2.2) 

The matrix Mij(s) consists of (n — l) 2 different functions, depending on the 
same set of arguments xi, . . . ,x n : 

Mij(s) = a 2 i-j[i -j + s], i, j = 1, . . . ,n - 1. (2.3) 



The main goal of this and next section is to transform the matrix ( |2.3j ) to 
such a form, that entries of a new matrix will be parameterized by a single 
function, depending on different sets of variables, (like V(xi,Xj) in (|2.1|) ) 

Mij^Dij, Dij = D(xi,Xj). (2.4) 

Here D(x,y) is a function of two arguments. The element depends on i 
and j only by means of its arguments and Xj. 

First, it is useful to rewrite the representation ( |2.2| ) in terms of a deter- 
minant of a matrix n x n . To do this, notice that Oin-j — 0, if j < n, so 
M n j = 5 n j[s] Y[m=i x m- Thus, we obtain 

n 

Qn{s) = Is]' 1 ]J x^detMjis), i,j = l,...,n. (2.5) 

m=l 

The r.h.s. of ( |2.5| ) is well defined for all s ^ and n ^ 0. These two cases 
should be considered separately. It is easy to see that for s = and ti / 
one have to understand ( |2.5| ) as a limit s —>■ 0, because the determinant is 
proportional to [s]. However, for n — 0, the original representation (|2.2|) is 
not well defined, while it is natural to define the determinant detoMjj = 1 
in ( |2.5|) . So, we obtain Qo(s) = [s]^ 1 for s / 0. On the other hand the 
case s = corresponds to the form factor of local field. In this case we have 
H = 0, and the form factor is equal to zero H o Q (0) = F = 0. Thus, we 
define Qo(s) = [s] _1 for s ^ 0. We do not define Qo(0), but we simply put 
F = for s = 0. 
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In order to study correlation functions we need to find the square of 
polynomials Q n (s). 

n 

Ql(s) = Is]' 2 II rt^nC jk , (2.6) 

m=l 

where 

n 

C jk = (M T ■ M) jk = Y / [i~J + s}[i-k + s]a2i-j(T2i-k- (2-7) 

i=l 

One can calculate the sum in ( |2.7| ) using an integral representation for ele- 
mentary symmetric polynomials 

^ = i/^n (z+a (2 - 8) 

Here the integral is taken in positive direction with respect to an arbitrary 
circle \z\ = p around the origin. Choosing the radius of the circle p > 1 and 
using ( |A.5| ) from Appendix A, we find 

n 

•^m ) 

r<., — I d 2 7 H!=l n-j+l n-k+l 

{ 2n+2+2s-j-k -2n-2-2s+j+k k-j 'j-k \ 
1 |_ 1 1 1 L _ (2.9) 
q 2 z\z 2 — 1 q~ 2 z\z 2 — 1 z\z\ — 1 z\z\ — 1 J 

In order to get a common denominator we make replacements of variables 
in the braces: Z\q x l 2 = W\, z 2 q 1 ^ 2 = w 2 in the first term; Z\q~ x l 2 = Wi, 
z 2 q~ l l 2 = w 2 in the second term; Z\q x l 2 = wi, z 2 q~ 1 ' 2 = w 2 in the third term 
and Z\q~ x l 2 = u>i, z 2 q l l 2 = w 2 in the fourth term. After simple algebra we 
arrive at 



where 



^ = (S5F> *» {2M) 



n 

n—t 



q q \ m=l m=l / 

(2.11) 
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The matrix C still depends on n 2 different functions Cj k . However, this 
matrix can be transformed to a more convenient form. Let us introduce 
matrix Aj k (it is studied in Appendix B) 



1 d n ~i n 



(2.12) 

i=0 



D jk = f d 2 w .n 2 / 2 2 2 —Y(w 1 ,x j )Y(w 2 ,x k ), (2.16) 



with a determinant 

n 

detA= Y[(x a -x b ). (2.13) 

a<b 

Instead of matrix C it will be convenient to introduce matrix D 

D = A T CA. (2.14) 
Determinants of matrices C and D are related by 

n 

det n C =U(x a - x b )- 2 det n D, (2.15) 

a>b 

The calculation of the explicit expression for matrix D in ( p. 14 ) reduces 
to the summation of the Taylor series (see (|B.8| )), so we have 

WiW 2 

where 

Y M =WL(* J£— Y (2.17) 

q — q \wq L/z + x wq L/z + x J 

and 

n 

J{w) = J] (m 1 ' 2 + x m )(wq- 1 / 2 + x m ). (2.18) 

m=l 

Taking the integral, for instance, with respect to w 2 (recall that |wiw 2 | > 1), 
we have after the symmetrization of the integrand 

D jk = ^/^( y ( w > x i) + Y(-w,x j ))(Y(w-\x k ) + Y(-w-\x k )). 

(2.19) 

Thus, we obtain a new representation for the square of the polynomial Q n (s): 

Ql(s) = • (2.20) 

[s] 2 n^n (x a - x b ) 2 

m=l a>b 

This brings us closer to (|2.1[). 
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3 Dual fields 

The entries of the matrix Djk are parameterized now by a single function 
D (|2.19 ). However, an element Djk, is still not a function of two arguments 



only, because of the product J(w) = rim=i( u "? 1/ ' 2 + ^m)^? _1 ' 2 + x m)- This 
product depends on all x m . In order to get rid of these products we introduce 
auxiliary Fock space and auxiliary quantum operators — dual fields. Dual 
fields are linear combinations of canonical Bose fields, see page 210 of . 
Let us define 

$iO) = qi(x) +p 2 (x), 

(3.1) 

where operators Pj(x) and qj(x) act in the canonical Bose Fock space in a 
following way 

((%,(£) = 0, Pj (x)\0) = 0. (3.2) 
Non-zero commutation relations are given by 

[Pi(x),qi(y)} = [P2(x),q 2 (y)} =£(x,y) = log ((a; + yq 1/2 )(x + yq~ 1/2 j). (3.3) 

Due to the symmetry of the function £(x,y) = £{y,x), all fields $j(x) com- 
mute with each other 

[%(x),$ fc (y)] = 0, j,k — 1,2. (3.4) 

However, despite of these simple commutation relations, the vacuum expec- 
tation value of the dual fields may be non-trivial, for example: 

(0|$i(x)$ 2 (i/)|0) = (0\p 2 (x)q 2 (y)\0) = £(x,y). (3.5) 

It is easy to show that an exponent of dual field acts like a shift operator. 
Namely, if / (<&i(y)) is a function of $i(y) then 



| J] e*^f($i(y)) |0) = (0| n e Mxm) f(qi(v)) |0) 

m=l m=l 

(0|/ U(y) + E ^ m ,y)) |0) = /(log J{y)). (3.6) 



m=l 
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Using this property of dual fields one can remove the products J(w) from 
the matrix Dj k . 
Let us define 

e *i(«0 / q s g-s \ 

Y(w, x) = — — — ZT75 - — , (3.7) 

q — q 1 ywq 1 '^ + x wq i ' z + x J 

and 

Djk = ^-f— (Y(w, x,) + Y(-w, x,)) (Y(w-\ x k ) + Y{-w-\ x k )) . (3.8) 
Then, due to (|3.6|), we have 

n 

det n D = (0| J] e* 2(xm) det n £>|0), (3.9) 



m=l 



det n D = (0|det„ (Dix^x^e^^ 2 ^) |0). (3.10) 



or 



The entries of the matrix Dj k depend on Xj and x k only, and they do not 
depend on other variables x m . Thus, we have presented the square of the 
polynomial Q n (s) in terms of a vacuum expectation value of a determinant 
of a matrix n x n, similar to one of the terms in the r.h.s of ( p.l|) . Entries of 
matrix D are parameterized by the single two- variable function D(x, y). Let 
us emphasize again that as an operator in the auxiliary Fock space D(x,y) 
belongs to an Abelian sub algebra. 



Besides the polynomial Q n (s) the form factor ( |1.4|) is proportional to a 
double product n">6 -^min (Pab) {%a + ^fe) -1 - in order to transform (fCTl) to 
(|2.1| ) it is necessary to factorize this product. To do this we introduce another 
dual field 

$oO) = q (x) +po(x). (3.11) 

As usual 

(0|g o (a;) = 0, p (x)|0)=0. (3.12) 

Operators qo(x) and po(y) commute with all pj and qj (j = 1,2). The only 
non-zero commutation relation is 

\po(x),q (y)]=Ti(x,y), (3.13) 
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where 



r)(x,y) = rj(y,x) = 2 log 



x 2 — y 2 



(3.14) 



Here we have used the fact that |F min (z) \ = \F min (—z) |. It is worth men- 
tioning also that the r.h.s. of ( |3.14j ) has no singularity at x = y, because 

F m m (x) has the first order zero at x = and F^ in (0) = (i sin (iir)' 
(see pi). Hence 



T](x, X) 



-2 log 



2 KB 

2x sin — F min (in) 



Newly introduced dual fields also commute 

[Q (x),$ (y)] = 0=[Q (x),$ j (y)]. 
However, due to the Campbell-Hausdorff formula, we have 



(oi n ^ 

m=l 

where 



|0) = U e ^ Xb) = \- n J] x m 2 ii 



a,b=l 



m=l 



i>b 



log — 

to Xj, 



2 / >-» 2 

x a x b 



A = 1 2 sin— F min 



27T 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



Combining the last formula and the representations fl2.20|) , (|3.10|) for Qn( s ) 
we find 



Qi(s) n 

a>b 



pp(0|det n f(x„x fc )|0). 



Here 
and 



X a -\- X\) 

V{xj,x k ) = D jk e^°^ + ^ ^ 
%{x) = %{x) + <5> 2 {x). 



(3.19) 



(3.20) 



(3.21) 



So we managed to represent a square of an absolute value of the form factor 
as a determinant, similar to one of the terms in the r.h.s of ( |2.1| ). In the next 
section we shall sum up all contributions of the form factors and obtain a 
determinant representation for a correlation function. 
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4 The determinant representation for a cor- 
relation function 

In the previous sections we have obtained the representation for a square of 
an absolute value of the form factor 

F n (/3 1 ,...,l3 n ) = (0\O(0,0)\Pi,---,Pn), (4- 1 ) 
in terms of a determinant 

\F n ((3 u ...,(3 n )\ 2 = j^\H n \ 2 \ n (0\det n V(x v x k )\0), (4.2) 
Here constants H n are equal to 

4 s i n eH 

Hin+i = H 1 ^, n , H 2n = H /j, n , n = — — -— , (4.3) 

and 

-L Ho = F = (0|O(0,0)|0), H 1 = F 1 = (0\O{0, 0)|/3). (4.4) 
[«] 

We have the representation for a correlation function of operators O in 
terms of form factors: 

oo . Jn a n 

(0|O(0,0)O(x,t)|0> = 5: / ZJ7^\U^---^n)?J[e^\ (4.5) 
where 

111 V 

e{x) = ^-(x + x- 1 ). (4.6) 
Substituting here (|4.2|) - (|43|) we arrive at a following representation 
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n=0 



o n! \ 27r 



xdet r 



!j^lf*l e -|W*jM«*)) 



\H \ 2 - \H^\^_ g d n s f |A/i|' 



n=0 



o ni 



2tt 



xdet r 



|0). (4.7) 



Both of these series have the form ( j2.1| ), so they can be summed up and 
written as determinants of integral operators (Fredholm determinants) 

(O|0(O,O)0(M)|O> = (Ol|^| |g ° |2 + l f l|2H " 1 det(/ + 7 £/) 



+ |if ° |2 ~ l 5 l|2N " 1 det(J- 7 ^)1 10), 



(4.8) 



where 



and 



(4.9) 



4 2 irB 
7 = — sm . 

' n 2 



(4.10) 

The determinant representation (|4.8|) is the main result of the paper, 
therefore we summarize here the basic definitions. 

The integral operators / ± act on a trial function f(x) as 



[(/ ± >yU)f\{x) = f(x) ± 7 / U(x, y)f(y)dy. 



The kernel U (x, y) is equal to 



U(x,y) 



D(x,y) 



(4.11) 



(4.12) 
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where 



D(x,y) 



and 



1 



dw 



w 



(Y(w, x)+Y(-w, x)) (Y(w-\ y)+Y(-w-\ y)) , (4.13) 



Y(w, x) 



3*1 (w) 



Q 



— q 1 I wg 1 / 2 + x u>g l / 2 + 



a; 



(4.14) 



The dual fields Qq{x) and $i(x) were defined in the section |3] (see (|3.1| ) 
and ( |3.11| )). The main property of these dual fields is, that they commute 
with each other, so the Fredholm determinants det(/±7?7) are well defined. 
Certainly det(I ±jU) are operators in auxiliary Fock space, but they belong 
to the Abelian sub algebra. On the other hand, the vacuum expectation 
value of these operators is non-trivial. It follows from commutation relations 
( |3.3| ), ( |3.13| ), that in order to calculate the vacuum expectation value, one 
should use the following prescription 



Mi 

lie 

a=l 



*o(z a ) 



M 2 

n 



3 *i(^) 



10) 



Mi Mi 

0=1 6=1 



Mi M 2 

> xi) n n e 

a=l 6=1 



Here 



and 



r}{x,y) = 2 log 



x 1 - 



€(x,y) = log((x + yg 1/2 )(x + yq 1/2 ; 



(4.15) 



(4.16) 



(4.17) 



Recall also, that the determinant representation (|4.8|) is valid for an ar- 
bitrary s. If s = 0, one should understand the r.h.s. as a limit s — > 0, taking 
into account that H = 0. 

Similar Fredholm determinant representations were useful not only in the 
theory of correlation functions ||-[|r2], but also in matrix models [P^]-[ J16|J . 
A work on determinant representation for correlation functions led to the 
discovery of a determinant formula for a partition function of the six-vertex 



model with domain wall boundary conditions [30]. In the paper [31] it was 



shown that this partition function satisfies Hirota equation. In the paper |32 



it was shown that the determinant formula for the partition function of the 
six-vertex model helps to solve a long-standing mathematical problem — to 
prove the alternating sign matrix conjecture. 
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5 Large r-asymptotic 



In this section we shall demonstrate, how one can find a long distance asymp- 
totic of a correlation function starting from the Fredholm determinant. We 
shall reproduce some known results. 

The kernel of the integral operator U(x, y) can be written in the form 



U(x,y) = j> dwP 1 (w,x)P 2 (w,y), (5.1) 
where projectors P are 

Pi(w,x) = —-(Y(w,x) +Y(-w,x))e 1 ^ o{x) - 1 2 e{x \ (5.2) 



ykvw\Jx 

P 2 (w, y) = -L(y(w~\ y) + Y{-w-\ y)) e ^°(3/H%). (5.3) 

Let us remind here that w integration goes along a large contour around 
zero in positive direction. A radius of the contour should be greater then 
1. The Fredholm determinants of the kernels of type ( |5.1| ) can be written as 
determinants of operators acting in the space of variables u w" 

det(I ± ~fU(x, y)) = det(J ± -fU(w u w 2 )), (5.4) 

where 



U(w 1 ,w 2 ) = / dxPi(wi, x)P 2 (w 2 , x) . (5.5) 
Jo 

The integral operator U(wx,w 2 ) acts on a trial function f(w) as 

[(I + U)f](wi) = f(wi) + <fiJ(w h w 2 )f(w 2 )dw 2 . (5.6) 



Consider the case r — > oo. Then the value of the integral in (|5.5|) can be 
estimated by means of a steepest descent method. The saddle point of the 
function 8(x) is x = xo = 1. Examination of commutation relations of dual 
fields ( |4.15| ) shows that dual fields can be considered as analytic functions in 
the vicinity of real axis. Hence, we can estimate the integral in ( |5.5|) as 



U(wi, w 2 ) = PxK, l)P 2 (w 2 , 1) I 4/ — + 0(r- 3 / 2 ) I . (5.7) 



mr 
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Thus, for the large r asymptotics the kernel U(wi,W2) becomes a one-di- 
mensional projector, and its Fredholm determinant is equal to 

det(/±7C7) -»• 1±7 £ dwU(w,w). (5.8) 

In order to calculate a vacuum expectation value of U(w, w) one can use 
prescription (|4.15 ), however it is better to write down the dual field $o(x) in 



terms of the original fields $o ( x ) — &o ( x ) + ^2 {%) • Then the contribution of 
the fields $o(l) gives 

(0|e* o(1) |0) = e H 1 > 1 ) = A" 1 . (5.9) 
To find a contribution of the fields $1(10) and $2(2^) we can use ( |3.6|) and 

(EH 



(0\U(w, w)\0) = ——J — fY^w, 1) + Y 1 (-w, 1) 
Sttzaw V rnr v 

x(y 1 (^- 1 ,i) + y 1 (-«;- 1 ,i)), (5.10) 



where 

{wq 1 ' 2 + l)(u;g- 1 / 2 + 1) / q s 



Y 1 (w,l) 



q — q 1 \ wq 1 / 2 + 1 wq' 1 / 2 + 1 



= [s] + [s - l/2]w. (5.11) 

After substituting this into ( |5.8|) it becomes clear that only a pole at w = 
contributes into the integral, so we arrive at 



dwU(w,w)\0) = [s] 2 \-\—e~ mr , (5.12) 



mr 



and hence 



det(J ± 7 LT) 10) -> 1 ± [s] 2 ? J— e^ mr . (5.13) 

A V Tnr 
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Finally, substituting this into (|4.8|) , and using explicit expressions for A 
( |3.18|) , fi Q4.3|) and 7 (|4 . 1 0|) we obtain the correct asymptotical expression 



(0|£>(0,0)£>(x,t)|0> -> 1^1! + \H 1 \ 2 (2Trmry l/2 e- mr . (5.14) 



Recall, that for the correlation function of local fields one should put Hq = 
= Hq/[s], therefore we see, that asymptotic formula ( |5.14fl is well defined 



for arbitrary s. If Hi = (for the stress-energy tensor), then (|5.14 ) gives 



a constant for an asymptotic. However, in this case one had to estimate 
the kernel U(wi,W2) more accurately. Namely, one should take into account 
explicit expression for corrections of order r~ 3 / 2 in ( |5.7|) . In this case the 
kernel U turns into two-dimensional projector and it is easy to show, that 
exponentially decreasing term behaves like exp(— 2mr). 



Summary 

We was able to sum up contributions of all the form factors and to obtain 
the closed expression for correlation functions ( |4.8| ). We believe that it will 
be possible to do not only for Toda models |28| but also for all models of 
integrable massive relativistic field theory. 
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A Elementary symmetric polynomials 

Consider an integral representation (|2.8|) for the elementary symmetric poly- 
nomials 

1 r dz n 
a^\x 1 ,...,x n ) = a k = — f Y[{z + x m ), (A.l) 

zmJ\z\= P z m=1 

where p is arbitrary positive. Notice, that the representation (|A.1|) holds for 
arbitrary integer k including k < and k > n. 

Here we are deriving auxiliary formula, which is used in the section ^| 
Namely, let us consider a sum 

n 

Tjk(a) = a2la 2i-j&2i-k, j, k = 1, . . . , n, (A. 2) 

i=l 

where a is an arbitrary complex number. One can extend the summation in 
l.21) from — oo to n. Then we have 



n 



Tjk{a) = J! a 2l <J2i-jV2i~k = ^Z^ 11 V) (52n-2l-j<y'2n-2l-k- (A. 3) 

i=~oo 1=0 



Using the integral representation (|A.1|) we find 



1 00 / d 2 za 2n ~ 21 n 

= J^nV ^ T U+j-n+121+k-n+l II (*1 + X m)(z2 + X m ). (A.4) 
l z/u J 1=0 J z l z 2 m=l 

We can choose the integration contour in such a way that \aziz 2 \ > 1 at the 
contour. Then one can sum up the series with respect to I: 

2n+2 ,, n-j+1 n ~k+l n 



ex."'" ' " p z z 

T ^ = 7^2- T d2z ^2^J i II (*i + X ™)(Z2 + x m ). (A.5) 
yzTYi) j a z 1 z 2 — i m=1 



The integrand contains only two simple poles aziz 2 = ±1, therefore one 
can take the integral with respect to z\ or z 2 and obtain a single integral 
expression for Tjk- 
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B Properties of the Vandermonde matrix 

Consider a Vandermonde matrix Wjk" 



W 



_fc-i 



jk — Zj , 



= l,...,n, 



with 



det n (W} fc ) = JJ(z a - z h ). 

a>b 



(B.l) 
(B.2) 



The inverse matrix W 1 can be written in a form 

1 d?- x n 



(j - 1)1 dxi 



7 n 

mj^k 



(B-3) 



i=0 



Indeed 



n n—1 ~,l jI n 

1=1 1=0 UJj m^k 



= n 

2=0 m ^ fc 



5jk- 



(B-4) 

Here we have used the fact, that the r.h.s. of (B^) is a Tailor series for the 
polynomial of the (n - 1) degree Um^k( x ~ z m){ z k ~ z^ 1 . 
In the section |] we used the matrix A, fc : 



1 d n ~^ n 
{n-j)\ dx n ' ^ k 



(B.5) 



x=0 



The determinant of this matrix is equal to 



det n A = J[(x a - x b ) det 



A? fc IT ( Xm ~~ X k) 
mj^k 



(B.6) 



It is easy to see that the matrix in the r.h.s. of ( |B.6| ) coincide with inverse 
Vandermonde matrix W^ 1 up to replacement x m = —z m and a permutation 
of rows. Thus, we obtain 



det n A = Y[{x a - x b ). 

a<b 



(B.7) 
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A calculation of products of the matrix A and matrices, containing powers 
of some complex numbers w, is simple. For example, deriving ( j2.16| ) we used 



g(^7y!^n. (x+ ^ 



x=0 



q-q 



s n n ( l/2\n-l Jn-l n 

* rn(».«- i/2 +^)E l! ^^^n(-+ 



m=l 



(n — 1)\ dx r 



„~s n n ( — l/2\n— I Jn—l n 



q-q 



m=l 



(n — l)\ dx r - 



7 



x=0 



x=0 



11 (wig 1/2 + x m )(w 1 q 1/2 + x m ) r 

m=l 



q-q 



-i 



Wiq 1 / 2 + Xj Wiq l l 2 + x 



j J 



• (B-8) 
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